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Optomechanical and electromechanical cavities have been widely used in quantum memories and 
quantum transducers. We theoretically investigate the robustness of opto(electro)-mechanical quan¬ 
tum memories against the noise of the control laser. By solving the Langevin equations and using 
the covariance matrix formalism in the presence of laser noise, the storing fidelity of Gaussian states 
is obtained. It is shown that, the destructive effect of phase noise is more significant in higher 
values of coupling laser amplitude and optomechanical coupling strength G. However, by further 
increasing coupling coefficient, the interaction time between photons and phonons decreases below 
the coherence time of laser frequency noise and the destructive effect of laser phase noise on the 
storing fidelity drops as well. 
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I. INTRODUCTION 


Optical quantum memories are one of the most impor¬ 
tant and essential equipments in quantum computers and 
quantum repeaters [l[. Several methods for preparation 
of these memories have been developed over the last few 
decades One of the most famous method is based 

on the Electromagnetically Induced Transparency (FIT) 
effect, which has been observed in lambda-shaped atoms 
il. Using this effect, by changing the coupling coeffi¬ 
cient between light and atom, the group velocity of light 
can approach zero which results in mapping quantum 
state of light to an ensemble of atoms, quantum state 
ini- The upper limit of the storage time is of the order 
of the decoherence time of the atom [HI- In recent years 
optomechanical systems have had many applications spe¬ 
cially as a macroscopic lab for observing quantum me¬ 
chanical effects. By adjusting the input laser detuning, 
the mechanical mode of these systems can be cooled close 
to the ground state [l^ - fTsj l . Recently, the FIT like effect 
has been observed in optomechanical cavities which is 
called optomechanically induced transparency (OMIT) 
[Tbl - fisjl . OMIT can provide a new approach to designing 
quantum delay lines and memories |l9l . [^ . Quantum 
state conversion between optical and mechanical modes 
of cavity has also been studied actively [2ll - [^ . It has 
been shown that, in the strong optomechanical coupling 
regime of operation, if the optomechanical coupling is 
quantum coherent, the optical and mechanical modes are 
exchanged with the frequency proportional to the op¬ 
tomechanical coupling coefficient [^, . This effect 

can be applied in the design of optical quantum memo¬ 
ries and transducers which transfer the quantum state of 
light from one frequency to the other . The stor¬ 

age time of optomechanical quantum memories is limited 
to i where 7 is the mechanical decay rate of the oscil¬ 
lator. Optomechanical memories and transducers have 
been studied theoretically as well as experimentally in 
many groups [H, [2§ - l^ . In reference [ 2 ^ , the storing 
fidelity has been obtained based on the quantum state 
diffusion method (QSD). As it is clear, the noise from 


the control (coupling) laser has a damaging effect on the 
storing fidelity. Fundamentally, the intensity and phase 
noises of the control laser are related by the uncertainty 
principle. In our analysis it is assumed that by stabiliz¬ 
ing the laser intensity, the intensity noise is negligible and 
phase noise becomes the main component of the control 
laser noise. The effects of laser phase noise in cooling 
rate of mechanical mode and entanglement, have been 
studied by many researchers [Sll - I^ . According to Rabl 
et al. if we consider the laser phase noise to be white 
noise, we won’t be able to cool the the mechanical res¬ 
onator to reach the ground state. However, experiments 
have shown the ground state cooling to be a possibil¬ 
ity. Due to the finite coherence time of laser frequency 
noise, the phase noise must be considered as a color noise 
with finite bandwidth png. The phase noise spectral 
density is equivalent to the noise spectral density of a 
low pass filter excited by white noise [s^. The cutoff 
frequency of the low pass filter is inversely proportional 
to the coherence time of the laser frequency noise. The 
flat spectral density, corresponds to the zero coherence 
time of input frequency noise (white noise). In this pa¬ 
per the effect of control laser phase noise on the storing 
fidelity of a Gaussian quantum state in the mechanical 
mode of the cavity, is investigated. It is shown that, 
the quantum memory fidelity destroy drastically in the 
presence of white phase noise. The effect of laser fre¬ 
quency noise coherence time on memory fidelity is also 
taken into consideration. In order to achieve the strong 
optomechanical coupling, the microwave electromechan¬ 
ical system that are easily cooled to temperatures below 
100 mK [ 3 ^ and an optomechanical Fabry- Perot system, 
are considered [i^. Effects of optomechanical coupling 
strength on memory fidelity in the presence and in the 
absence of laser phase noise, is studied. It is shown that 
by raising the intracavity amplitude, the destructive ef¬ 
fect of phase noise increases. This paper is organized as 
follows: in section II, the theoretical model of the system 
is explained and the dynamics of the variables consid¬ 
ering the effect of phase noise, are obtained. The noise 
model is also introduced in this section. The method of 
obtaining the storing fidelity of a Gaussian state and the 


2 


storage protocol is described in Section III. Section IV is 
devoted to the numerical results and analysis of the effect 
of phase noise on the storing fidelity of electromechani¬ 
cal and optomechanical systems, and finally the paper is 
concluded in Section V. 


II. THEORETICAL MODEL 

We consider an optomechanical system for storing the 
quantum state of light. This system can be electromag¬ 
netic cavities such as Fabry-Perot [i^, microsphere or 
microtoroid cavities with one mechanical mode [ 2 ^, or 
electromechanical systems based on microwave resonance 
circuit [ 3 ^ . The schematics of a Fabry-Perot cavity with 
one movable mirror and a microwave circuit are shown in 
Figs. 1(a) and 1(b), respectively. The cavity is driven by 
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Figure 1: (color online) Schematics of the following systems: 
a) An optomechanical Fabry-perot cavity consisting one vi¬ 
brating mirror. The cavity is driven by a laser from the left 
mirror, b) A microwave cavity. 

an intense coupling laser beam with frequency ujl , which 
is red detuned from the cavity resonance frequency ujc- 
The radiation pressure force of the driving field leads to 
displacement of the movable mirror, changing the res¬ 
onance frequency of the cavity and coupling the intra¬ 
cavity field to the mechanical mode of the mirror. In 
microwave cavity, the displacement of the movable plate 
of capacitor due to the radiation pressure force of the mi¬ 
crowave field, changes the capacitance C and resonance 
frequency of the microwave circuit (wq = where L 

is the inductance of the microwave resonator [ 3 ^ . High 
quality factor microwave resonators can be reached by 
employing super conductors. In the resolved sideband 
limit (when the mechanical resonance frequency LOm ex¬ 
ceeds the cavity decay rate k), if the coupling laser is 
tuned to the lower optomechanical sideband, in the ro¬ 
tating wave approximation (RWA), the linearised interac¬ 
tion Hamiltonian between mechanical and optical mode 
of the cavity is Hint = —hG{aH + am) which is beam¬ 
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splitter-like Hamiltonian [ 2 ^ [ 2 ^ [ 2 ^. a and b are the 
linearised bosonic annihilation operators for cavity opti¬ 
cal photon and mechanical phonons, respectively. G is 
the coupling coefficient between the optical and mechan¬ 
ical modes, which is proportional to the driving field. In 
the absence of any optical and mechanical dissipation, 
a{t) and b{t) can be derived easily: 

b{t) = b(0)cos{Gt) + ia{0)sin{Gt) 

a(t) = a(0)cos{Gt) + ib{0)sin(Gt) (1) 

As it can be seen from Eq. o, the quantum state of 
the mechanical and optical modes can transform into one 
another (up to a phase factor) with period This 

conversion is the basis of the design of quantum memory 
by opto(electro)-mechanical systems. Initially we should 
prepare the optical mode of the cavity in the state |' 0 ) we 
wish to store it in. By using the optomechanical coupling 
f pulse {At = ^), the optical state of interest is mapped 
to the mechanical mode. After the storage time r is up, 
the second optomechanical coupling ^ pulse will map the 
state that is stored in the mechanical mode to optical 
mode and we can retrieve it [1^ [ 2 ^. In practice due 
to the existence of optical and mechanical energy decay 
rates and input noises, the initial pure state converts to 
mixed state and the fidelity of the final state decreases. 
In order to have a quantum memory, the optomechanical 
coupling must be larger than both the optical decay rate 
K and the thermal decoherence rate of the system 'jth = 
NmJ {ith, K < G), where and 7 are the mean thermal 
phonon number of reservoir and mechanical damping rate 
of the oscillator, respectively. If this condition is satisfied 
the optomechanical coupling is quantum coherent [ 2 ^. 
Our study is concentrated on the effect of the driving 
laser phase noise on the retrieved state fidelity in the 
quantum coherent regime of operation. 

The Hamiltonian of the system in a frame rotating at 
the driving laser frequency can be written as: 

H = HuJmbH + HAqo,^ a — hg^a^a(b + w ) 

+ fhEL(t)(e*‘"(‘)a^-e-*‘"(‘)a) (2) 


where Aq = luc — wl is the detuning of coupling laser 
from the cavity resonance frequency. (p{t) is the fluctuat¬ 
ing phase noise of the driving laser. EL{t) = EL+SEL{t) 
is the complex amplitude of the driving laser inside the 
cavity which is composed of one deterministic term [E^] 
and one fluctuating term {SEL{t)) which arises from fluc¬ 
tuations in the amplitude of the laser. The deterministic 


term is related to the input power P from E^ 



In a stabilized laser source, the amplitude noise of the 
driving laser is negligible compared to the phase noise, 
so we neg lect the amplitude noise of the laser in this 
paper [^. go = —^^Xzpf is the single photon optome¬ 
chanical coupling strength where L is the length of the 

cavity and Xzpf = 2 mui zero-point fluctuation 

amplitude of the mechanical oscillator, with m being the 
effective mass of the mechanical mode. 
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A. Quantum Langevin equations 


Quantum Langevin equations of the system can be de¬ 
rived by adding the corresponding noise and damping 
terms to the Heisenberg equations: 


= -iAod + igod(b + b'^)+£JLe"'^'^*^ - ^a + \/Kdin 


db 

dt 


igoa'a- -b- 




( 3 ) 


where din a-nd are input vacuum and thermal noise 
operators with two times correlation functions given by 


mm-. 


(a*„(t)dj„(t)) = {l + Nc)S{t-t) 
{d\n{t)d,n{t)) = Nc6{t-t) 

{b^n{t)blni^)) = {l + Nm)S{t-t) 
{Hn(t%nit)) = Nm6(t-t) (4) 

Here A^m(C) = mean thermal 

e _1 

phonon(photon) number of the thermal(optical) reservoir 
in the frequencies uJmi^c) at temperature T and Kb is 
Boltzmann’s constant. At optical frequencies Nc — 0. 
The statistical properties of phase noise will be deter¬ 
mined later. Using the transformation d —?► we 

will be in the rotating frame with instantaneous fluctu¬ 
ating frequency ip{t) |36l |. 


dd 

dt 

db 

dt 


—iAgd — iipd + igadib + + El — ^d + y/Hdin 

-iuJmb + igod'd “ 2 ^ (5) 


The standard method is used for linearising the non¬ 
linear equations ([5]). Operators are written as the sum of 
steady state values and fluctuating terms with zero mean 
value {d ^ as + Sd, b ^ f3 + Sb). The steady state values 
(fixed points) are obtained by setting the time deriva¬ 
tives of operators to zero {as = ,(Ao-go^+/ 3 U)+f ’ ^ ^ 
2 

i^°+i )• The algebraic equation for fixed point as is a 
polynomial equation of degree 3, which has 2 stable and 
1 unstable root [H, . In order to obtain a higher cou¬ 

pling coefficient, the highest intensity with a stable fixed 
point is chosen. The input laser phase is opted such that 
as is real. The dynamics of the fluctuating operators 
around these steady states are as follows: 


= —iASd — iipas + igoa{Sb + Sb^) 

— '^dd ~\- y/Hdin 
dSh 

— = -iuJrnSb + igoas{6d + 6d^) --^Sb + y/jbini6) 


where A = Aq — go{l3 + f3*) is the effective cavity laser 
detuning. We work in the resolved sideband regime and 


have a good cavity {k <C ujm)- In the rotating wave ap¬ 
proximation (RWA), by choosing A = uim-, the dynamical 
equations will be written as: 


dd 

dt 

db 

dt 


—iuimd — iipas + iGb — —d + y/Hdin 
-ioJmb + iGd --b + y/jbin (7) 


where G = goas is the coupling coefficient between opti¬ 
cal and mechanical field and d, b will be used as the lin¬ 
earised operators from now on. We define quadratures of 

each mode as: Xi = (-^- 5 ^),Pi = ) where i = 1 , 2 

and Ai = 6 , A 2 = d. Using Eq. (O one can obtain the 
evolution of the quadratures in the matrix form: 


— = QX + 8 

dt 

where X G 3?^ is the vector of the quadratures : X = 
[xi,f)i,X2:P2Y ^nd Xin G is the input noise opera¬ 
tors vector: X„ = A/KP 2 ,m - 

agifiY' . XLin,Pi,in{i = 1,2) are defined by fi,i„ = 
) with Bi = bin,B 2 = din and 

Q is given by: 


Q 


( 

2 



_1 

2 

-G 

0 


0 

G 

K 

”2 




0 



(9) 


B. The noise model 

The statistical properties of Xi^in,Pi,in{i = 1,2) are 
defined by using Eq. (SD- The phase noise of the laser 
is responsible for the Laser linewidth T l- By taking the 
statistical properties of the laser phase noise as a white 
noise, the correlation function of the frequency noise will 
be {(p{t)(p{t)) = 2 r L5{t — t) and a flat noise spectrum is 
obtained S'i^(w) = 2Tl- However due to the finite non 
zero correlation time of frequency noise, the frequency 
noise spectral density can not be flat. In other words, 
it has a finite bandwidth color noise. In the simplest 
model, the spectrum is equivalent to a low pass filtered 
white noise with the following spectrum and correlation 
function [3l|: 


impit)) = rL7ce-"^l‘-‘l (10) 

-L “r —T 

7c 

Here 7 ^ is the inverse of frequency noise correlation time 
or cutoff frequency in the laser frequency noise spectrum, 
so that the frequency noise is suppressed at frequencies 
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a; > 7c- By increasing 7c, the correlation time decreases 
and the frequency noise starts reaching the white noise. 
This frequency spectrum corresponds to the differential 
equation (f(t) + = £(t) where e(t) is a Gaussian 

random variable with white noise correlation function 


(e(t)e(t)) =2'Y^rLS{t-t) 


( 11 ) 


In the following calculations, it is more convenient to 
use input white noise rather than color noise, so we 
define an auxiliary operator tl) = ip and add its cor¬ 
responding differential equation (tpit) + lc‘P{t) = ^( 0 ) 
to our four coupled differential equations of quadra¬ 
tures. The quadratures vectors X and Xin are ex¬ 
tended to Xe = [x\,Pi,X2,P2t'^V £ ^11*1 Xin,E = 

[V7ii.m, in, v^i2.m,\/KP2,*n,e]’^ G 5 ?’", respec¬ 
tively and Eq. dS]) is changed to = Qe^e + 
with: 


Qe = 


In the following section, the method to hnding the fidelity 
of the final state with the initial state, in the presence of 
coupling laser is explained. 


III. GAUSSIAN STATE FIDELITY 


so if pf is the rotated version of pi, the storing fidelity 
will be 1. 

Two useful parameters fih = 2 y/det{Vi + Vf) — 1 and 

= {{^^) - (I/)) • - (if)) is introduced 

to make the fidelity formula more straightforward [2^: 


F = 


I 


l + Uh 




(14) 


The parameter h/j which is the effective number of ther¬ 
mal quantas, originates from input noises and heats the 
state throughout the process. The parameter A shows 
the effect of optical and mechanical damping rates on 
the fidelity of the state. In order to study the fidelity 
of the retrieved state, the expectation value of optical 
quadratures and also covariance matrix at initial and fi¬ 
nal states, must be calculated. The initial covariance 
matrix is obtained easily by knowing the initial state of 
the optical and mechanical mode. The evolution of the 
covariance matrix can be derived from Eq. (I8|) (see Ap- 


/ 

7 

2 


0 

-G 

0 ^ 

pendix A): 


Cu’m 

7 

2 

G 

0 

0 



0 

-G 

K 

2 

COrn 

0 

(12) 


G 

0 


K 

2 



V 

0 

0 

0 

0 

-7c / 

where N is 


dV ^ 

— QeV + VQe 


N 


(15) 


color noise as (See Appendix A): 


N = 


/ ^ 0 0 0 

' 0 I 0 0 

0 0 f 0 

0 0 0 “ 


(16) 


V 0 0 0 0 272 rL / 


For the sake of simplicity, it is assumed that the initial 
state of the cavity is a pure Gaussian state. At this step 
we are going to find the fidelity of the optical state at 
time t with a given initial state, while the optomechanical 
interaction is switched on. Our approach is the same as 
the one used in reference [ 1 ^ . The fidelity between initial 
and final state is given by F = Tr{piPf), where pi and 
Pf are the density matrices of the initial and final states, 
respectively. In Wigner representation, the fidelity can 
be written as E = tt J^^d^Wi{^)Wf{^), where ^ G 

is the vector of the optical quadratures ^ = [x2,P2]’^- 
The Wigner function of Gaussian states can be written 
in terms of covariance matrix of quadratures V, where 
the symmetrized covariance matrix is defined by Vij = 
^{XiXj + XjXi). Due to the bilinear form of interaction 
Hamiltonian, the final state remains Gaussian. So after 
some mathematical manipulation, a relation for Gaussian 
state fidelity is achieved: 


F = 


+ Vf 




(13) 


Indices i and / stands for initial and final states, respec¬ 
tively. It is noted that the fidelity must be maximized by 
employing an orthogonal transformation in phase space. 


7 and R are defined by 7 = 7(1 -I- 2Nm) and R = 
k(I -I- 2 Nc), respectively. Eq. (ITSl) is an inhomogeneous 
differential matrix equation that can be written as a sys¬ 
tem of 25 coupled ordinary differential equations. Due to 
the symmetry in the covariance matrix, only 15 of these 
equations are independent. By solving these equations 
by ordinary methods like RungeKutta, the covariance 
matrix of the final state can be derived. However, the 
analytic solution of this equation is as follows: 


V{t) = e^®V( 0 )e' 5 «‘ 



It is assumed that the average value of input noises are 
zero. Averaging over Eq. ([ 81 ), the input noises are elimi¬ 
nated and the equation is reduced to : = Q{X(t)), 

which has the following solution: 

(i(t)) = eO‘(A(0)) (18) 


Note that due to elimination of noise, we don’t have 

to use the extended version of X and Q. Using Eqs. 
(HZl),® and m , the optical state fidelity at time t 
with initial state can be derived. As it has been men¬ 
tioned before, the terminal fidelity should be maximized. 
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Using the unitary transformation to rotate the pa¬ 

rameters on the phase space, the fidelity is maximized, 
where Qr is given by: 


Qr 


( 0 

0 

0 0 
0 0 
Vo 0 


0 0 0 \ 

0 0 0 

0 0 
0 0 

0 0 0 / 


(19) 


If the initial optical state is not a Gaussian state, such 
as fock state and any finite linear combination of fock 
states, it would be difficult to find analytic expression 
for the fidelity. Generally the initial density matrix of 
the system can be written as p(0) = Po ® Pm , where 
indices o and m indicate optical and mechanical modes, 
respectively. The evolution of density matrix of the sys¬ 
tem can be obtained by the numerical solution of the 
Lindblad master equation in the presence of o ptic al and 
mechanical dissipations and input noises 1^. The 
final density matrix can be used in F = Tr{p{0)p(t)) to 
give the fidelity between initial and final states. 


A. Storing Protocol 


The storing fidelity depends on the input state. Mem¬ 
ory fidelity can be defined as the average of fideli¬ 
ties over the possible input states. In this paper our 
study is concentrated on the squeezed coherent state: 
|a,r) = D{a)S{r)\0) where D{a) = “ and 

S{r) = are displacement and squeezing op¬ 

erators, respectively. It is assumed that r is the real 
squeezing parameter. By setting r = 0, we will have a co¬ 
herent state. The initial state of the oscillator is assumed 
to be a ground state |0), which means that the oscilla¬ 
tor has been cooled. This assumption is almost possible 
by using a different cooling method (l2l - [l5l| or using an 
optomechanical coupling ^ pulse before the storing pro¬ 
tocol, which is actuall y c ooling the oscillator near the 
quantum ground state [23 ■ If the cooling process hasn’t 
been done, the initial state of the mechanical mode would 
be the thermal state. The storing protocol is composed 
of three steps: 

1 ) At initial time, the average of the quadratures are 
(A(0)) = [0,0, Re{a), Im{a)]^. The initial covariance 
matrix is: 


U(0) 


/ 1 0 0 0 0 \ 

0 10 0 0 

0 0 e-^’- 0 0 

0 0 0 0 

V 0 0 0 0 1 / 


( 20 ) 


If the initial state of the oscillator hasn’t been cooled, 
two elements of the covariance matrix V (Un and V 22 ), 
are changed to 1 -I- 2Nm, where Nm is the mean ther¬ 
mal phonon number as defined before. By switching on 
the coupling laser for the duration time G = ^ (write 


pulse), the state of the optical mode is almost swapped 
with the mechanical mode. The quadratures average and 
covariance matrix can be obtained from Eqs. (11711 and 

m- 

2 ) The coupling laser is switched off and the optical 
state is stored in the mechanical mode for a time r. In 
this step, the Q matrix is changed to Q', where it is the 
same as matrix Q except that the elements containing G, 
go to zero. 

3) After the storage time is up, by switching on the 
coupling I pulse again (read pulse), the state that is 
stored in the mechanical mode, is mapped to the optical 
mode and can be retrieved from the cavity. The second ^ 
pulse must be tt radians out of phase relative to the first 
one to cancel the accumulated phase during the process 
[ 2 ^. The schematic figure of cooling, writing and reading 
pulses are shown in Fig. 2. 



Figure 2: (color online) A schematic outline of coupling ^ 
pulse for quantum memory. The first pulse applies for cool¬ 
ing the mirror near ground state. This pulse will change the 
state of initial optical mode (vacuum state) with thermal me¬ 
chanical mode. After a while, the intracavity thermal optical 
mode transfers to the reservoir and the input Gaussioan state 
enters the cavity at the beginning of the write pulse. After 
the storage time r, the reading ^ pulse with an opposite sign 
of the writing pulse, enters the cavity. 

The final quadratures average and covariance matrix 
is obtained from: 


(X/) = ( 21 ) 

and 

^Qe'’’ qQe 

E^^e'’’ ]\[^Q E E 

Jo 

+ f (22) 

Jo 

where matrix Q" and Q'^ is obtained by replacing G by 
—G in matrix Q and Qe, respectively. (Due to the phase 
difference between read and write pulses). The fidelity 
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is obtained by Eq. m- Effects of the environment are 
considered in cavity and mechanical decay rates and also 
in thermal phonon Nm and photon Nq numbers. Effects 
of control laser noise on the memory fidelity are given by 
presence of laser bandwidth El and frequency spectral 
density bandwidth jc- 

IV. NUMERICAL RESULTS 

To investigate the effect of coupling (control) signal 
noise on the memory fidelity, numerical calculations are 
done for a microwave cavity and an optomechanical cav¬ 
ity in the following subsections. 

A. microwave cavity 

To consider the effect of coupling signal phase noise on 
the storing fidelity in an electromechanical system, the 
parameters are chosen as those employed in references 
[2^ 1^. These parameters correspond to a microwave 
cavity. In the next subsection the optomechanical Fabry- 
Perot cavity, is considered. The resonance frequency of 
the mechanical mode is ^ = 10.69 MHz and the cav¬ 
ity decay rate is chosen to be ^ = 170 kHz [l^. We 
are working in the strong coherent coupling regime of 
operation <C k <C G) and so we choose the largest 
value of electromechanical coupling constant G that has 
been achieved in experiments to date (G = 0.05 x uj^n = 
27r X 534.5 kHz) [^. The initial value of electromag¬ 
netic mode that we are interested in storing, is the co¬ 
herent state I a) with a = 1. The single photon coupling 
strength go is 27r x 230 Hz which is obtained by opti¬ 
mizing both and Xzpf to achieve large values of go 
[ 3 ^. The mean thermal phonon number is taken to be 
Nm = 3, which corresponds to the temperature 1.7 mK 
and the storage time is set to r = 64 = 0.95 gs. 

By changing the mechanical quality factor of the system 
Qm = ^ (wm is fixed while 7 is changed), the storing 
fidelity E, nt and A, are calculated and plotted for differ¬ 
ent values of microwave control linewidth T^ and cutoff 
frequencies 7 c in Figs. 3(a), 3(b) and 3(c), respectively. 
By increasing the noise intensity in the electromechani¬ 
cal system, the parameters Uh and A are increasing which 
leads to the drop in fidelity. The amplitude and band¬ 
width of phase noise spectral density are proportional to 
the control linewidth T^ and cutoff frequency 7 c, respec¬ 
tively. Hence, the fidelity decreases by increasing T^ and 
7 c. As it is expected, Fig. 3(a) shows that, by increas¬ 
ing the cutoff frequency, the fidelity will mainly degrade. 
However, for the small values of the phase noise param¬ 
eters like El = 1 kHz and 7 c = 100 kHz (dash-dot 
green line in Fig. 3(a)), the phase noise will not affect 
the fidelity significantly. The limit 7 c —>■ 00 corresponds 
to the white noise, which shows very low fidelity nearly 
zero (purple line with marker o). Fidelity is an increas¬ 
ing function of mechanical quality factor Qm- Due to 
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Figure 3: (color online) Variation of a) Storing fidelity, b) 
heating parameter fih ans c) damping parameter A, versus 
mechanical quality factor Qm = Qim = 2t:x 10.69 MHz), 
for En = 0 (solid red line), En = 1 kHz, 7c = 100 kHz (dash- 
dot green line), Ei, = 1 kHz,^^ = 0.5 MHz (dashed blue line) 
and Ei, = 1 kHz,'yc = 10 MHz ( 7 c ^ 00 or white noise) 
(purple line with marker o ) (white noise hasn’t been shown 
in (b,c)). The input state is a coherent state |a) with a = 1, 
the thermal phonon number is Nm = 3, and the storage time 
is set to r = 

UJrn 


decoupling of mechanical oscilattor from its environment 
in high values of Qm, the effective number of thermal 
quantas that heat the state decreases, which results in 
decreasing the parameter fih (see Fig. 3(b)). Further¬ 
more, by decreasing mechanical decay rate 7 , the damp¬ 
ing parameter A decreases (see Fig. 3(c)), and as a result, 
the fidelity increases. 

For Qm > 10^, the storing fidelity for low values of 
noise parameters (dash-dot green line in Fig. (3.a)) is 
more than F > 0.8, which is suitable for experimental 
investigations. However, for low values of Qm, where the 
mechanical decay rate becomes large, the final state will 
thermalize faster and the fidelity declines. 

For a given mechanical quality factor of the oscillator 
Qm = 360000 (which corresponds to 7 = 186.57 s“^), the 
effect of the signal control linewidth E ^ and noise band¬ 
width 7 c in the fidelity is calculated and presented in Fig. 
(4.a) and Fig. (4.b), respectively. As it is expected, the 
fidelity is a decreasing function of signal linewidth E l , 
while in terms of cutoff frequency 7 c, it decreases and 
approaches saturation corresponding to the white noise 
value. Variation of fidelity as a function of cutoff fre¬ 
quency 7 c, for three values of E^ = 1 kHz, 5 kHz, 10 kHz 
, is presented in Fig. (4.b). 

The control phase noise is present only during the write 
and read process, while the environmental noises of op¬ 
tomechanical system exist in all the read, write and stor¬ 
age times. Fig. 5 shows the storing fidelity as a function 
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Figure 4: (color online) a) Variation of fidelity as a function of 
control linewidth Fl- b) The variation of fidelity as a function 
of the control cutoff frequency 7 c for Fi = 1 kHz (solid blue 
line), Fi, = 5 kHz (dash-dot green line) and F^ = 10 kHz 
(dashed red line). The input state is a coherent state |a?) with 
a = 1. = 360000. 

of the storage time r, for two values of reservoir temper¬ 
ature in the presence and in the absence of control phase 
noise effect. As the storage time increases, the duration 
time that mechanical mode interacts with thermal noise 
increases, which results in reducing fidelity. Reservoir 



Figure 5: (color online) Storing fidelity as a function of 
the storing time r for different values of reservoir temper¬ 
ature, laser linewidth F l and cutoff frequency yc '■ T — 
1.7 mK,TL ~ 0 (solid red line), T = 1.7 mK,TL ~ 
1 kHz, 7 c = 300 kHz (doted blue line), T = 0.01 mK, Fl = 0 
(dashed green line) and T = 0.01 mK,TL = 1 kHz,yc = 
300 kHz (dash-dot black line). The input state is a coherent 
state I a) with a = 1. 

temperature plays an important role in the storing fi¬ 


delity. If the reservoir temperature is 1.7 toA', the mem¬ 
ory will store the initial coherent state up to 0.4 ms with 
fidelity more than 0.8 (in the absence of control noise ef¬ 
fect) and with hdelity more than 0.7 (in the presence of 
control noise effect). By increasing the reservoir temper¬ 
ature to 0.01 K, the storage time drops to 0.05 ms for 
an appropriate fidelity. As it is expected. Fig. 5 shows 
that by increasing the storage time, the effect of control 
phase noise on hdelity decreases relative to other envi¬ 
ronmental noises and the curves considering the control 
phase noise effect, reach the corresponding curve with¬ 
out this effect. That is because the control phase noise 
inhuences memory only during the read and write pulses 
and doesn’t have any effect on the hdelity in the storage 
time. 

To consider the effect of coupling coefficient G = aggo 
on the hdelity of memory, an electromechanical system 
with the following parameters is taken into consideration. 
The mechanical quality factor of the oscillator is hxed 
and is equal to Qm = 360000. By changing the coupling 
coefficient G from 0.02 x uJm to 0.05 x uim, the hdelity of 
storing the coherent state |a) for storage time t = 
is calculated and represented in Fig. 6(a). The coupling 
coefficient G is adjusted by as- So, Fig. 6(a) also shows 
variation of memory hdelity versus as- As already ex¬ 
pected, in the absence of control phase noise (solid red 
line in Figs. 6(a) and 6(b)), the storing hdelity increases 
by increasing the coupling coefficient between microwave 
and mechanical modes. In this situation, the interac¬ 
tion time tg decreases as G increases, which results in 
less thermal noise affecting the system in the interaction 
time. However in the presence of control phase noise, the 
damaging effect of phase noise is more effective in higher 
coupling coefficient G, which results in the decrease of 
the storing hdelity in high values of G. This phenomenon 
originates from the coefficient ag of the phase noise term 
in Eq. ©• So, by increasing the input coupling (con¬ 
trol) intensity, despite of the drop in interaction time tg, 
the ehect of input control phase noise on the system gets 
stronger, which prevents the hdelity from increasing. Ac¬ 
cording to the mean value theorem, there is a maximum 
in the hdelity versus ag in the presence of control phase 
noise, which can be observed obviously in Fig. 6(a). 

Let us now increase the coupling coefficient G more 
than has been achieved in experiments. Fig. 6(b) shows 
the variation of memory hdelity versus coupling coeffi¬ 
cient, where G is changed from 0.02 x uJm to ujm- It is 
easy to show that in the strong coupling regime of oper¬ 
ation, an optomechanical system is equivalent to a tune 
hlter for the control phase noise, and its center frequency 
can be determined by the coupling coefficient G. Hence 
the control phase noise can be modelled by cascading a 
low pass hlter with cutoff frequency 7 c and a band pass 
hlter with a central frequency determined by the coupling 
coefficient G (see Fig. 6(c)). Increasing ag increases the 
noise intensity and consequently decreases the memory 
hdelity. By increasing ag even more, the tune hlter cen¬ 
tral frequency moves away from the cutoh frequency of 













G/k 


G/k 



(c) 


Figure 6 : (color online) a, b) Storing fidelity as a function 
of ^ for two ranges of variation of G (a,b) , in the ab¬ 
sence and in the presence of laser phase noise effect: F = 0 
(solid red line), F = 1 kHz,^c = 100 kHz (dashed green 
line), r = 1 kHZj'jc = 200 kHz (dash-dot blue line) and 
F = 1 kHz,'yc = 300 kHz (dotted black line), c) Schematic 
description of cascading a low pass filter with cutoff frequency 
7 c and a band pass filter with center frequency determined 
by the coupling coefficient G for two values of coupling coef¬ 
ficient Gi cc cTsi and G 2 oc as 2 - Ois 2 > Qsi. The storage time 
isr= 


low pass filter and the noise intensity decreases due to 
normal mode splitting. In other words, the read and 
write duration time decreases below the coherence time 
of control frequency noise and as a result the memory 
fidelity increases. The minimum observed in Fig. 6 (b) is 
related to this effect. For large Og, the effect of control 
phase noise is negligible and fidelity is affected by the 
optomechanical system’s environmental interactions. 

In order to compare the storing fidelity of the squeezed 
state with the coherent state, the squeezed state fidelity 
|Q:,r-) for four values of r = 0,0.2,0.5,0.8, as a function 
of Qm is calculated and shown in Fig. 7(a). The storage 
time is set to r = 0.95 and other parameters are 
chosen like those are used in Fig. 3. The values of F^, 7 c 
is taken F/, = 1 kHz and 7 c = 10 kHz, respectively. 
Fig. 7(b), shows the variation of fidelity as a function of 
squeezing parameter r, when mechanical quality factor 
is fixed to Qm = 360000. It is clear that the fidelity 
decreases as parameter r increases while the initial state 
becomes more and more non classical. By growing the 
squeezing parameter r, the initial state would be more 
sensitive to input thermal and optical noises [ 2 ^. 



Figure 7: (color online) a) Storing fidelity as a function of 
mechanical quality factor Qm with F= 1 kHz, 7 ^ = 10 kHz 
for different values of squeezing parameter r: r = 0 (solid 
red line), r = 0.2 (dash-dot green line), r = 0.5 (dashed 
blue line) and r = 0.8 (dotted black line), b) Variation of 
fidelity as a function of squeezing parameter r for fixed vane 
of Qm = 360000. The storage time is set to r = 0.95 /rs. 


B. Optomechanical cavity 

The optomechanical cavity that we consider is a Fabry- 
Perot cavity with strong coupling between optical and 
mechanical mode [40| . The cavity and mechanical de¬ 
cay rates are ^ = 215 kHz and ^ = 140 Hz, re¬ 
spectively. The mechanical frequency of micromirror is 
= 947 kHz and the single photon coupling strength 
is 5 o = 27r X 1.91 iJz. The input coupling Nd:YAG 
laser with power 11 mW and wavelength 1064 nm enters 
the cavity and makes optomechanical coupling coefficient 
G = 27rx 229.81 kHz. In the absence of laser phase noise, 
the hdelity of storing the coherent state |q;) with a = 1, 
for storage time 0.95 /rs, with thermal phonon number 
Nm = 3, will be obtained 0.789. By considering the effect 
of phase noise as a laser linewidth F^ = 1 kHz and cutoff 
frequency 7 c = 1 kHz, the storing fidelity decreases to 
0.28 which is an undesirable fidelity. In order to compare 
the fidelity of microwave and optomechanical systems, 
the storing fidelity for some values of noise parameters 
as a function of the mechanical quality factor (by fixing 
ujm and changing 7) is presented in Fig. 8. The storage 
time and thermal phonon number have been chosen to 
be similar to Fig. 3 (r = 0.95 , Nm = 3 ). 

Comparing Fig. 3 with Fig. 8, we hnd out that the 
degrading effect of phase noise on the fidelity of optome¬ 
chanical system is more signihcant than that in the mi¬ 
crowave system. This could be explained by considering 
the fact that, the single photon coupling strength 50 of 
our microwave system (27r x 230 Hz) is more than that 
one in optomechanical system (27r x 1.91 Hz). So in 
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noise as a control line width 1 kHz and cutoff frequency 
7 c = 300 kHz, the storage time with appropriate fidelity 
decreases to 0.05 ms. By increasing the input control 
intensity, the optomechanical coupling G increases and 
the interaction time between photons and phonons de¬ 
creases, which results in memory fidelity gain, while at 
the same time due to the increase in intracavity field am¬ 
plitude, the effect of phase noise elevates which prevents 
the fidelity from increasing. By comparing the electrome¬ 
chanical and optomechanical systems in the range of ex¬ 
perimentally achievable parameters [H, , it was found 

that, the system with higher single photon coupling go, 
is more robust against phase noise degradation effects. 


Acknowledgments 


Figure 8: (color online) Storing fidelity as a function of me¬ 
chanical quality factor Qm, for Fn = 0 (solid red line), 
Fn = 1 kHz,^c = 0.5 kHz (dashed green line), Fn = 
1 kHz,'yc = 1 kHz (dash-dot blue line). Storage time and 
thermal phonon number are r = 0.95 fis and Nm = 3, respec¬ 
tively. 

order to achieve the same value of coupling coefficient 
G, the lower intracavity amplitude as is needed in mi¬ 
crowave system. The reduction in as leads to the re¬ 
duction of phase noise effect in microwave system. On 
the other hand, the coupling coefficient G in microwave 
system (G = 27r x 534.5 kHz) is more than that in op¬ 
tomechanical system (27r x 229.81 kHz) which gives the 
shorter interaction time between light and matter in mi¬ 
crowave system and hence the fewer amount of thermal 
and phase noise, will affect the system. 


V. CONCLUSION 

The effect of laser phase noise on the fidelity of elec¬ 
tromechanical and optomechanical quantum memories 
have been studied. The phase noise of the laser behaves 
like a color noise which is included in the dynamics of 
the equations by means of an axillary variable with in¬ 
put white noise. It is shown that the damaging effect 
of laser phase noise on memory fidelity is proportional 
to the intracavity field amplitude as- For the electrome¬ 
chanical system that we have considered in the strong 
coupling regime of operation, the storing fidelity of the 
coherent state |a) for an interval r ^ 1/xs, drops from 
0.95, in the absence of control phase noise effect, to 0.66, 
in the presence of control phase noise, with a linewidth 
of Tl = 1 kHz and a cutuff frequency of jc = 0.5 MHz. 
By decreasing the cut off frequency to 7 c = 10 kHz, the 
memory fidelity increases significantly and the degrad¬ 
ing effect of noise disappears. It has been shown that 
by minimizing the noise effect, the coherent state can be 
stored up to 0.4 ms with a fidelity more than 0.8, in me¬ 
chanical mode. While by considering the effect of phase 
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Appendix A 

The dynamic equations of quadratures vector are ob¬ 
tained from the extended version of Eq. (|S|). Here we 

define variable: Vi = {XX'^). It is easy to show that 
the covariance matrix is related to Vi through equation 
V = . Using Eq. ([5]), the following equations are 

derived: 

XX^ = QeXX^ + X,„X^ 

+ (Ai) 

which gives us the relations: 

QeVi + Uigi + 

+ QV/’ + (i„i^ + (A2) 

By adding these two equations, the dynamic equation 
of the covariance matrix is achieved: 

dv 

dt + 2 

(A3) 

The last term of Eq. (IA3I) is the diffusion matrix N. Solv¬ 
ing Eq. ®, the quadratures can be written according to 
the input noises as: 

X{t) = e'3‘A(0) + [ e'5(*-‘)A„(/) (A4) 

Then the diffusion matrix will be given as: 

N = i(A + A^) (A5) 

where 

nt 

A= dt[D{t,t) + D{t,t)'^]e^^^*-^'> (A6) 

Jo 


dVi 

dt 

dV^ 

dt 
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where D{t,t) = {Xin{t)Xj^{t)) and it is derived easily 
using the correlation functions of the input noises: 


diffusion matrix N is obtained by employing Eqs. 
(IA61) and (IA7|) : 


D{t,t) = 


5{t — t) 


/ 

7 

17 

0 

0 



/ 

7 

A 

0 

0 

0 

0 \ 


—ij 

7 

0 

0 
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0 

1 

0 

0 

0 


0 

0 

K 

in 

0 

II 

fT 

< 


0 

4 

0 

K 

0 

0 

V 

0 

0 

—in 

K 

0 



0 

0 

4 

0 

K 

0 

0 

0 

0 

0 

87c"ri y 


V 

0 

0 

0 

4 

0 

j 


(A 8 ) 


where 7 = 7(1 + 2Nm) an d k = nil + 2Nn). The 
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